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ABSTRACT 

We explore the onset of star formation in the early Universe, exploiting the observations of high-redshift 
Lyman-break galaxies (LBGs) and Lya emitters (LAEs), in the framework of the galaxy formation scenario 
elaborated by Granato et al. (2004) already successfully tested against the wealth of data on later evolutionary 
stages. Complementing the model with a simple, physically plausible, recipe for the evolution of dust atten- 
uation in metal poor galaxies we reproduce the luminosity functions (LFs) of LBGs and of LAEs at different 
redshifts. This recipe yields a much faster increase with galactic age of attenuation in more massive galaxies, 
endowed with higher star formation rates. These objects have therefore shorter lifetimes in the LAE and LBG 
phases, and are more easily detected in the dusty submillimeter bright (SMB) phase. The short UV bright 
lifetimes of massive objects strongly mitigate the effect of the fast increase of the massive halo density with 
decreasing redshift, thus accounting for the weaker evolution of the LBG LF, compared to that of the halo mass 
function, and the even weaker evolution between z~6 and z w 3 of the LAE LF. The much lower fraction 
of LBGs hosting detectable nuclear activity, compared to SMB galaxies, comes out naturally from the evolu- 
tionary sequence yielded by the model, which features the coevolution of galaxies and active nuclei. In this 
framework LAEs are on the average expected to be younger, with lower stellar masses, more compact, and 
associated to less massive halos than LBGs. Finally, we show that the intergalactic medium can be completely 
reionized at redshift z w 6 - 7 by massive stars shining in protogalactic spheroids with halo masses from a few 
10'°Mq to a few 10' ' Mq, showing up as faint LBGs with magnitude in the range -17 < M1350 < -20, without 
resorting to any special stellar initial mass function. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: high-redshift - galaxies: luminosity 
function - intergalactic medium 



1. INTRODUCTION 

The impressive recent advances in the observational studies of high redshift galaxies and quasars have brought us up to the 
outskirts of 'dark ages', when the amount of neutral intergalactic hydrogen starts becoming significant and tentative indications 
of the presence of metal-free population III stars begin to show up. This may give us some confidence that we are approaching 
the onset of the galaxy formation history and provides key information on physical processes driving the early galaxy evolution 
and on their influence on the intergalactic medium (IGM). 

Different techniques and selection criteria have been exploited to find high-z galaxies. A most efficient method is the Lyman 
break dropout technique (Steidel & Hamilton 1993), first applied to select z « 3 galaxies (Steidel et al. 1996; Lowenthal et al. 
1997) and subsequently extended to redshifts of up to « 10 (Steidel et al. 1999, 2003; Dickinson et al. 2004; Giavalisco et al. 
2004; Ouchi et al. 2004; Bunker et al. 2004; Vanzella et al. 2005, 2006; Bouwens et al. 2005; Thompson et al. 2006; Bouwens 
& Illingworth 2006a,b). 

Very deep narrow band imaging at selected redshift windows (Cowie & Hu 1998) has also been successfully used to search for 
Lya emission up to z « 7 (Ouchi et al. 2003; Kashikawa et al. 2006b; Shimasaku et al. 2006; lye et al. 2006). 

While these techniques are sensitive to unobscured or moderately obscured objects, very dusty distant galaxies with extremely 
high star formation rates (SFRs) have been detected by SCUBA 850/im surveys (see e.g. Chapman et al. 2005; Aretxaga et al. 
2007). Substantial numbers of passively evolving high-z galaxies have also been found (Fontana et al. 2006; Grazian et al. 2007). 

The different selections emphasize different aspects of the early galaxy evolution that need to be integrated in a unified view. In 
this paper we present a simple physical model featuring an evolutionary link between Lya selected (LAE), Lyman break (LBG), 
submillimeter bright (SMB), and passively evolving galaxies, establishing quantitative relations between luminosities and basic 
physical quantities such as stellar and halo masses, star formation rates, galactic ages, metal abundances and dust attenuation'. 
The model extends to earlier evolutionary phases the approach by Granato et al. (2004), according to which the evolution of 
massive galaxies is tightly linked to the growth of the supermassive BHs at their centers. The model has successfully accounted 
for an extremely large amount of observational data on galaxies and quasars (Granato et al. 2004, 2006; Cirasuolo et al. 2005; 
Silva et al. 2005; Lapi et al. 2006; see Appendix A for details). 

An extension of the model to the earliest evolutionary phases, when metals were first produced, requires a specific treatment 
of dust attenuation, which plays an important role in the derivation of high-z galaxy properties, since even small amounts of dust 

' We follow the convention, increasingly adopted in the recent literature, of using the term 'attenuation' rather than 'extinction' to indicate the difference in 
magnitudes between the observed and intrinsic starlight, which is the final result of the interplay between the dust properties and the complex relative geometry 
of stars and dust. The terra 'extinction' is reserved to the idealized situation in which a single source is attenuated by a thin screen of dust. 
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Fig. 1 . — Left panel: SFR as a function of the galaxy age for different halo masses Mh = lO'^M© (solid line), IO"M0 {dot-dashed line), lO'^M© {dashed 
line), and lO'^ Mq {dotted line), virialized at z = 6. Note the step-like cutoff in the SFR for the two higher masses, due to the quasar feedback (see Appendix A 
for details). Right panel: same for the stellar mass. 



Strongly affect the rest-frame UV emission. For example, a correction for dust attenuation corresponding to E{B-V) ~ 0.2 would 
increase the estimate of the SFR at z « 6 by a factor of about 10. The presence of significant attenuation in high-z LBGs has 
indeed been reported by many authors (e.g., Steidel et al. 1999; Shapley et al. 2001, 2006; Huang et al. 2005; Ando et al. 2006; 
Yan et al. 2006; Burgarella et al. 2006; Eyles et al. 2007; Stark et al. 2007b), not surprisingly since the dust formation in early 
galaxies is expected to occur on very short timescales (Morgan & Edmunds 2003). 

The outline of the paper is the following. In § 2 we present an overview of the model and show that, with a simple recipe 
for the evolution of dust attenuation during the earliest phases of the galaxy lifetime, accounting for the luminosity vs. E{B-V) 
correlation found by Shapley et al. (2001), it nicely reproduces the observational estimates of the LBG luminosity functions 
(LFs) at different redshifts. In § 3 we investigate the absorption effects on the Lyman continuum emission, the evolution of the 
Lya LF, and the contributions of young galaxies to the cosmic re-ionization. Finally, in § 4 we discuss our results and summarize 
our conclusions. 

Throughout the paper we adopt a flat cosmology with matter density Q,m = 0.3, Hubble constant Ho = 70 km s~' Mpc~', and 
normalization of the mass variance (j% = 0.8. Unless otherwise specified, magnitudes are in the AB system. 

2. OVERVIEW OF THE MODEL 

The rest-frame UV emission is proportional to the SFR M*. This is in turn determined by a number of physical processes 
(such as cloud-cloud collisions, cooling, fragmentation, energy inputs from quasars, supernovae [SNae] and stellar winds, etc. 
etc.) which occur on different linear and time scales. Even with the more sophisticated numerical techniques, one has to resort 
to many simplifying assumptions to derive the star formation history, since these processes are extremely complex and occur on 
scales well below the currently achievable resolution of A^-body+hydro simulations (Iliev et al. 2006). 

In this paper we exploit the physical recipe elaborated by Granato et al. (2004), which relates the average SFR to the distribution 
of gas and dark matter (DM) in proto-galaxies, taking into account the effects of the energy fed back to the intra-galactic gas 
by SN explosions and by accretion onto the nuclear, supermassive black hole (BH). The model envisages that during or soon 
after the formation of the host DM halo, the baryons falling into the newly created potential well are shock-heated to the virial 
temperature. The hot gas is (moderately) clumpy and cools fast especially in the denser central regions, where rapid mergers 
of a fraction less than 30% of the total mass yield strong bursts of star formation that can reach up to thousands solar masses 
per yr in massive halos at high redshift. Such a huge star formation activity also promotes the storage of the cooled gas into a 
reservoir around the central supermassive BH, eventually leading to accretion onto it. The ensuing SN explosions and the nuclear 
activity feed energy back to the gaseous baryons, and regulate the SFR and the BH growth. These mutual energy feedbacks 
actually reverse the formation sequence of the stellar component of galaxies compared to that of DM halos: the star formation 
and the buildup of central BHs are completed more rapidly in the more massive haloes, thus accounting for the phenomenon now 
commonly referred to as downsizing. The basic equations of the model and their analytical solutions are reported in Appendix A. 
The model predicts the time evolution of the SFR, of the metal abundance Z (taking into account the chemical enrichment of the 
primordial infalling gas and the gas outflow due to energy injection by SNae and active nuclei), of the stellar mass (see Fig.[T]l, 
and of the mass stored in the central supermassive BH. 

When the star formation and the chemical evolution history of a galactic halo of mass Mh have been computed, the spectral 
energy distribution (SED) as function of time from extreme UV to radio frequencies is estimated through the GRASIL code 
(Silva et al. 1998). Coupling these results with the halo formation rate, we can obtain the LFs of galaxies as a function of cosmic 
time. At z > 2, the halo formation rate can be reasonably well approximated by the positive part of the derivative of the halo mass 
function with respect to cosmic time (e.g., Haehnelt & Rees 1993; Sasaki 1994). The halo mass function derived from numerical 
simulations (Jenkins et al. 2001) is well fit by the Sheth & Tormen (1999, 2002) formula, that improves over the original Press & 
Schechter (1974) expression, which is well known to under-predict by a large factor the massive halo abundance at high redshift. 
Adopting the Sheth & Tormen (1999) mass function Nst(Mh,T), T denoting the cosmic time, the formation rate of DM halos is 



Dust in primeval galaxies 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



Oo oo * 



o 




0.0 0.1 0.2 0.3 

E(B-V) 



0.4 



Fig. 2. — Correlation of the intrinsic (attenuation-corrected) rest-frame absolute UV magnitude MJjjq with the color excess E(B—V). The open circles are 
the data by Shapley et al. (2001) at z « 3, and the dashed line shows the average observational relationship. The black dots show the expectations of our 
model using the attenuation eq. (3), for halo masses in the range IO'^Mq ^M// ^ 4 x 10'^ Mq sampled in intervals A logMu = 0.2, and for ages in the range 



2 X 10 yr < ? < A/bmst that correspond to > 70% of the star-forming periods. The relation E(B-V) ' 
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o/ll by Calzetti et al. (2000) has been adopted. 



Nst{MhJ), (1) 



where a = O.lQl and p = 0.3 are constants obtained by comparison with A^-body simulations, <t(Mh) is the mass variance of the 
primordial perturbation field, computed from the Bardeen et al. (1986) power spectrum with correction for baryons (Sugiyama 
1995) and normalized to erg w 0.8 on a scale of 8/i~' Mpc, and 5^(0 is the critical threshold for collapse, extrapolated from the 
linear perturbation theory. 

Briefly, the model requires 10 parameters (aU of which are constrained within rather small ranges by independent data or 
physical arguments) to describe the physical processes ruling the star formation, the BH growth and their feedback. With the 
parameter values reported in Table 1 of Lapi et al. (2006), the model fits the available data on the LFs of protospheroidal galaxies 
at different redshifts and in different spectral bands, and the local relationship between central BH mass and spheroidal luminosity 
(hundreds of data points). Among others, the model reproduces the 850 pm extragalactic counts (Coppin et al. 2006; Scott et al. 
2006) with their redshift distributions (Chapman et al. 2005; Aretxaga et al. 2007), and correctly predicted (Silva et al. 2004, 
2005) the fraction of high-z galaxies in deep 24 /im surveys (Perez-Gonzalez et al. 2005; Caputi et al. 2006; see Magliocchetti et 
al. 2007), quantities that all proved to be very challenging for competing models. 

The evolving LFs of optical and X-ray selected quasars at z > 1.5 are also well fitted by adding 3 parameters, namely the 
visibility times in the optical, Afvis = 5 x 10^ yr, and in X-rays, Afvis = 3 x 10^ yr, and a dispersion of 0.3 dex around the mean 
relationship between the halo and the central BH mass derived from the model (see Lapi et al. 2006). 

A summary of the model results and of their comparisons with observational data is presented in Table 2 of Lapi et al. (2006). 
To sum up, the model traces very well the formation and coevolution of quasars and of their spheroidal hosts. 

2.1. The UV luminosity function ofLBGs 
The absolute magnitude at A « 1350 A, including attenuation, can be written as 

£-1350 ^ . , 



Mi35o = 51.59-2.5 log- 



ergs 



Hz" 



-2.5 log 



Mq yr- 



+A 



1350 



(2) 



where L1350 ~ 1 .2 x 10^^ erg s ' Hz ' is the monochromatic luminosity for = \Mq yr ' (Kennicutt 1998, with the calibration 
converted to the Romano et al. (2002) initial mass function [IMF] adopted here, see Appendix A). 
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Fig. 3. — Attenuation (left panel) and attenuated magnitude (right panel) at 1350 A as a function of galactic age, computed using eq. (5) and the quantities 
M*(t) and Z{t) provided by our model. The shaded area shows the range probed by the observed LFs. Lines are for different halo masses virialized at j = 6, as in 
Fig.E] 

The dust attenuation is expected to be proportional both to some power of star-forming gas mass (or equivalently of the SFR, 
see eq. IIA5II ) that reflects the dust column density, and to some power of the metallicity, that reflects the fraction of metals locked 
into dust grains (see also Jonsson et al. 2006). As shown by Fig.|2] a simple, average relation between UV attenuation, SFR and 

consistent with expectations from the model, provides a very good fit of the luminosity-reddening relation found by Shapley et 
al. (2001). Variations of the normalization coefficient by 30%, of the exponent of the SFR by ±0.15, and of the exponent of the 
metallicity by ±0.2 are still compatible with the data. 

Using eq. (|3]l and the quantities M*(r) and Z(f ) provided by our model (see Appendix A for handy approximations), we compute 
the evolution with galactic age of the attenuation for LBGs with different halo masses, see Fig. |3] (left panel). The attenuation 
is predicted to be small for LBGs hosted by low-mass halos (My < lO^'M©), while it quickly increases with time for higher 
masses, so that more massive galaxies are UV-bright for shorter times. The more massive LBGs spend approximately 90% of 
their burst time in an interstellar medium optically thick to their UV emission. The right panel of Fig. [3] shows the predicted 
evolution of UV luminosity with galactic age, including the attenuation given by eq. (O, at z = 6 for several halo masses. 

The UV LF of LBGs can be estimated coupling the UV luminosity and halo mass derived from the our model (cf. Fig. O with 
the halo formation rate. The rest-frame UV LF at a cosmic time T is then 

^350,7^)= / dMH—^— r[Mi35o|MH,r] , (4) 

J AMndT dMi350 

where T[Mi35o|M//, T] is the time spent at magnitudes brighter than M1350. 

As illustrated in Fig. |4] good agreement with the observational data is obtained. Neglecting attenuation would increase the 
prediction at the bright end by more than an order of magnitude. In fact dust attenuation is drastically shortening the time 
T[Mi35o|M//, r] [eq. dUl] for large halos (cf. Fig.O. In this context it is worth noticing that the unattenuated LF can be depressed 
also by assuming short star formation timescales (around 80Myr at z 6), as advocated by Stark et al. (2007a). On the other 
hand, such short timescales would imply stellar masses and stellar ages much lower than observed; Eyles et al. (2005, 2007) 
find stellar masses around 1 - 3 x 1O'**M0 and ages of hundreds Myrs. These findings are quite naturally reproduced in our 
framework. 

The approximate proportionaUty of the initial (galaxy age f = 0) SFR with the halo initial gas mass which, in turn, is proportional 
to the halo mass [eq. ( IA6H implies that the highest luminosity tail of the LF is associated to the most massive galaxies. These 
are however very rare at high-z and, furthermore, their UV-bright phase is very short because of the quick increase of their 
attenuation (Fig. |3]l; they are thus easily missed by the available surveys, covering very small areas. On the other hand, low mass 
galaxies are too faint to be detectable. Hence, according to the present model, the available data only sample the halo mass range 
1O'"M0 <Mh< IO'^Mq, with SFRs of « 1 - lOOMgyr"', and stellar masses ranging from « 10^ to few x IO'^Mq. The 
strong dependence of A 1350(f) on halo mass has clearly a key role in shaping the UV LF. 

The information on the LBG space densities at z > 7 is scanty. Fig.|4]shows that we expect a drastic decrease from z = 7 (dashed 
line) to z = 10 (dotted line). Bouwens & Illingworth (2006b) estimate that there are 1 -4 sources brighter than M^go « -19.7 at 
7 < z < 8 on an area of 5.8 arcmin^; we predict 2 sources. Bouwens et al. (2005) obtained an upper limit of 3 galaxies at z ~ 10 
on an area of 15 arcmin'^, to a magnitude limit of 28.5 in //160; our model yields 1. 

3. THE EFFECT OF DUST ON THE LYMAN-CONTINUUM OF LBGs AND ON LAEs 



3.1. Escape fraction of ionizing photons 
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Fig. 4. — Model luminosityfunctions of LBGs at 1350 A for z = 3 (dot-dashed line), z = 5 (solid line), z = 7 (dashed line), and z = 10 (dotted line), based on 
the attenuation given by eq. }3). Data at z ~ 3 — 4 are from Steidel et al. (2001; empty circles) and Yoshida et al. (2006; empty squares); at z ~ 4-5 are from 
Iwata et al. (2001 , filled circles) and Yoshida et al. (lOOd, filled squares); at z ~ 5 — 7 from Bouwens et al. (2005; stars). The Steidel data have been converted 
from 1700 A to 1350 A setting M1350 ~ M1700 + 0.2, as suggested by Bouwens et al. (2005) on the basis of the composite spectrum of z ~ 3 LBGs. 

The relative importance of dust attenuation and HI absorption in determining the escape fraction of ionizing photons and the 
Lya luminosity depends on the complex geometry of the distribution of gas, dust and stars. We assume that most of the ionizing 
and of the Lya photons originate in the very central regions of the galaxy and see the same amount of dust as the photons at 
the reference wavelength of 1350 A; we extrapolate the 1350 A attenuation to shorter wavelengths adopting the attenuation law 
proposed by Calzetti et al. (2000). This yields A 1350 « IIE(B-V), Ame ~ l.lAogo, and, extrapolating from 1200 A to 912 A 
(cf. Draine 2003), A912 « I.6A1350. 

The observed luminosity at 912 A is given by 

robs_rint ^912^912^912 /^x 
^912 -^912 /hi /dust/iGM ' y-') 

where Lg"2 is the intrinsic luminosity at 912 A, f^\^ and /^^Jj^^j ~ ^-^912/1.08 jj^g fractions of 912 A photons surviving HI and 
dust absorption within the galaxy respectively, and /j^^ is fraction surviving absorption in the IGM along the line of sight. The 
product X /j^p « /esc is the fraction of 912 A photons emerging at the galaxy boundary. 
Following Inoue et al. (2006) we define the ratio R^sc between the Lyman continuum and UV luminosity at the galaxy boundary 

as 

T obs 

_ ^912 /--s, 
"esc - ^912 ■ \»J 

^1350 iiGM 

The factor fiQ^, which in principle should appear in this equation, is very close to unity and is therefore usually omitted. 

Rest-frame ultraviolet spectroscopic observations in the Lyman continuum region for samples of 29 and 14 LBGs at z w 3 
have been obtained by Steidel et al. (2001) and by Shapley et al. (2006). While Steidel et al. (2001) found R^,^ w 0.22 ±0.05, 
Shapley et al. (2006) found an average value about 4.5 times lower, /?esc ~ 0.05, yet with a substantial variance and with 2 objects 
showing significant emission below the Lyman limit, corresponding to /?esc ~ 0.35 and 0.22. Note that these authors determine 
/iGM b^sed on a sample of quasar spectra at the same average redshift of the LBGs, and/or performing numerical simulations of 
the line-of-sight absorption. A large scatter of the observed ratio is expected, since the UV and the ionizing radiation are emerging 
from central regions of galaxies with quite complex gas and dust distributions. In fact, regions of strong star formation in local 
galaxies exhibit line of sights with quite variable dust attenuation and Lya emission. In the following we adopt /?esc ~ 0.15 as a 
reference value, but we explore the range 0.1 < R^sc ^ 0.2. 

Considering that Ljjg,, = ^050 /Jus?' (recall that /^f" « //g^" w 1) and setting R^t = Lj," 2/^050, the escape fraction writes 

r _ ^esc ^1350 

/esc — „ /dust ■ V ' / 
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Fig. 5. — Top: escape fraction of ionizing photons as a function of galactic age. Middle: fraction of ionizing photons absorbed by HI within the galaxy as a 
function of galaxy age. Bottom: fraction of photons escaping the galaxy at the Lya wavelength. The lines correspond to different halo masses virialized at z = 6, 
as in Fig.[T] In all panels we have adopted the attenuation eq. (3), the Calzetti attenuation law and Rcsc = 0.15. The bumps in the solid lines are just numerical 
artifacts. 

The ratio between the intrinsic luminosities at 912 A and at 1350 A for continuous star formation depends only on the IMF and 
on the metal content. For the Romano et al. (2002) IMF adopted here (see Appendix A for details), we find Rmt ~ 0.3 (see also 
Leitherer et al. 1999; Steidel et al. 2001; Inoue et al. 2005). In Fig. |5] (top panel) we plot /esc as a function of galactic age, 
computed according to the above equation and assuming f^^^^ ~ 1350/1 -08. jj. apparent that for massive halos (Mh ^ lO'^M©) 
the escape fraction is negligible over most of the burst duration. 

On the other hand, the low mass (i.e. low intrinsic luminosity) LBGs have large escape fractions and can thus release substantial 
amounts of ionizing photons (Fig.|6] left panel). The emission rate of ionizing 912 A photons at the galaxy boundary reads 

^912= ^ «5.4x 10«/e. ( tAt) (8) 
hvgn \Mq yr ' / 

The contribution of these galaxies to the cosmic reionization is estimated in Appendix B. The filling factor of HII regions is 
computed as a function of redshift (see Fig.|6] right panel) for our reference attenuation [eq. ©I and Re^c ~ 0.15. It turns out 
that LBGs are able to completely reionize the universe since z « 7. The corresponding optical depth to electron scattering is 
Tes « 0.07, within the range of the estimate from WMAP three year data (res = 0.09 ±0.03, Page et al. 2007). Note that lowering 
7?esc to 0.1 and 0.05 would move the completion of reionization to z « 6.6 and 5.8, respectively. 

Most of the ionizing photons are produced by galaxies with halos masses in the range 10'°'^ M0 < Mh S lO^ '^M© at z < 8. 
The contribution of massive galaxies is depressed by their low number density at the relevant redshifts and by their internal dust 
absorption. On the other hand, the contribution of the numerous low mass galaxies is damped by the decrease, driven by stellar 
feedback, of their SFRs. LBGs in the above mass range have very low attenuation and contribute to the faint part of the UV LF 
(-20 < M1350 < -18) and to the bright end of the observed Lya LF at z = 6.5 (see below). 

3.2. Lya emitters 

The Lya line typically arises from ionizing photons absorbed by nearby hydrogen gas. If the fraction of absorbed ionizing 
photons within the galaxy is (1 -/hJ") and 2/3 of them are converted into Lya photons (Osterbrock 1989; see also Santos 2004), 
one gets the unattenuated Lya luminosity: 

iLyo=^^12(l-./S')/^^Lya. (9) 
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Fig. 6. — Left: evolution of the ionizing luminosity at 912 A as a function of the galactic age; lines correspond to different halo masses virialized at z = 6, as in 
Fig.[T] Right: evolution of the filling factor of ionized region in the Universe with the redshift. In both panels we have adopted the attenuation eq. ^3), the Calzetti 
et al. (2000) attenuation law and ficsc = 0.15. 

The interstellar dust attenuates Lj^^^ by a factor f^l^', furthermore, only a fraction f]^Q^ of Lya photons survives the passage 
through the IGM. Therefore the luminosity we see is 

L?/„ « 5.6 X 10« [j^) (1 -/S^)/]u^s?/iG^M erg s"' , (10) 

for the adopted IMF; for a Salpeter IMF the normalization should be reduced by a factor of about 1 .6. 
Since the fraction of ionizing photons surviving HI absorption is /j^p k, /esc//dusr eq. (|7]l yields 

p fl350 
j912 _ ^esc 7 dust 

/hi - p ^912 • ^^^> 
"int ./dust 

The fraction /j{J^ depends on the attenuation at 1350 A and on the attenuation law used to extrapolate it to 912 A. The middle 

and bottom panels of Fig. |5] show the absorbed fraction (1 -/hJ^) and its product by f^l^ as a function of the galactic age for 
the attenuation eq. Q and the Calzetti et al. (2000) attenuation law; note that for the dust surviving fraction we have assumed 
/<kist — e"'*^/' It is apparent that LBGs in large halos, Mh > lO'^M©, are LAEs only for a short time (less than 3 x 10^ 
yr). The figure also illustrates the difficulty of inferring the SFR from the Lya luminosity: even for relatively low mass halos 
{Mii < 10" Mq) the effect of HI absorption and dust attenuation may be very large. This implies that the Lya emission from 
LBGs and in general the statistics of LAEs at high redshifts are important absorption/attenuation probes. 

Exploiting the relationship between the Lya luminosity, the SFR and the absorption [eq. (fTOl)l. the LFs of LAEs at high 
redshifts are estimated in the same way as the LEG LFs, and are given by an equation strictly analogous to eq. (|4|i. Our results 
are presented in Fig.|7](left panel). 

It is clear from the above that the LF estimates rely on some quantities that are not very well constrained observationally. Still, 
it may be worth noticing that the data on the LF at z « 4.9 and z « 5.7 are consistent with a transmission factor fl^^ ~ 1, for 

7?esc ^0.15, while the LF at z « 6.4 seems to require f^^^ ~ 0.5. This may suggest that the ionization state of the IGM is changing 
between z ~ 5.7 and z ~ 6.4 (see also Kashikawa et al. 2006b). We need, however, to keep in mind that the determination of 
fj^^ is quite complex and implies assumptions on the relative velocity of the emitting galaxies and the IGM, on the presence of 
galactic winds, of neutral hydrogen, etc. (see Haiman 2002, Santos 2004 and Dijkstra, Lidz & Wyithe 2007 for comprehensive 
discussions). On the other hand, a fit of the observed LF with f^\^ « and ignoring dust attenuation yields a conservative lower 
hmit of > 0.1 -0.4 at z w 6.5 (see Dijkstra, Wyithe & Haiman 2007). 

Figure Q also shows our predictions for the LAE LF at z ~ 8, when the Lya line lies in the / band (see Barton et al. 2006), 
assuming fl^^ « 0.5. However, if indeed the reionization is due to LBGs, we expect that the IGM ionization level can be 
significantly lower at this redshift, and the dimming of the Lya luminosity substantially stronger. 

A clear-cut prediction of our model is that the Lya luminosity of the most massive galaxies decreases abruptly after few 10^ 
yr, implying that, with rare exceptions, bright LBGs should exhibit relatively low Lya luminosity (see Fig. |7] right panel). A 
deficiency of objects with large Lya equivalent width among bright LBGs has recently been reported by Ando et al. (2006). 

According to the model, most of the LAEs with Lya luminosities sampled by the available surveys (2 x 10^^ erg s"' < LLya S 
3 x lO"^-' erg s"' at z > 5 - 6) are hosted by galactic halos within a rather narrow mass range 5 x lO'^M© < Mh ^ 5 x lO^M© 
and shining only for a short time < 8 x 10^ yr Fig. [T]shows that the corresponding stellar masses range from a few 10^ M© to 
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Fig. 7. — Left: Model Lya LF for z = 4.9 {dashed line), z = 5.7 {dot-dashed line), z = 6.4 {solid line), and z = 8.2 {dotted line). Observational upper and lower 
limits at redshift z = 6.4 are from Kashikawa et al. (2006b; filled triangles), data at z = 5.7 are from Shimasaku et al. (2006; empty circles) and from Murayama 
et al. (2007; empty squares), and data at z = 4.9 are from Ouchi et al. (2003; empty stars). We have adopted the attenuation eq. (3)> 'he Calzetti et al. (2000) 
attenuation law, ifesc = 0.15, /j^jJJ = 1 at both z = 4.9 and z = 5.7, and /j^j^ = 0.5 at both z = 6.4 and z = 8.2. Right: hya luminosity as a function of the galactic 

age for Rcsc = 0.15 and /I'q^ = 0.5. The shaded area shows the range probed by the observed LFs. The Hues correspond to different halo masses viiialized at 
z = 6, as in Fig.[T] 

a few 10^ M0, in nice agreement with the findings by Filkenstein et al. (2007) and Pirzkal et al. (2007). Indeed, Gawiser et al. 
(2006) noted that 80% of the objects in a sample of LAEs at z « 3.1 have the right colors to be selected as LBGs. They also 
inferred an average SFR around 6Mq yr"' and an average mass in stars M,t « 5 x IQ^Mq, assuming a Salpeter IMF. Applying 
the correction by a factor of about 2 to account for the different IMFs, our model yields SFRs and stellar masses consistent with 
these estimates. Gawiser et al. (2006) also claim that the dust attenuation is on the average minimal (Ay < 0.1) in their stacked 
sample. According to our model (cf. Figs.[T]|3] and|7]i, the attenuation is almost negligible for low stellar mass objects, since both 
the SFR and the metal abundance are small. Only for objects with M^, > 3 x 10*^ and age larger than 5 x 10^ yr the attenuation 
becomes non negligible, A1350 > 0.3. The dust in these objects may play an important role in determining the EW of the Lya 
(Filkenstein et al. 2007). However the effect is very sensitive to the spatial distribution of dust, stars and HI. 

Lai et al. (2007) found larger masses, from a few xlO^M© to a few x10'"Mq, for 3 LAEs at z « 6. The corresponding 
ages vary from several Myr to several hundreds Myr, depending on chemical composition. It has to be noticed that these authors 
analyzed the 3 objects with the largest fluxes at 3.6 and 4.5 ^m. As a consequence we expect that the selection is biased in 
favor of the objects with the highest stellar mass. In our framework this selection corresponds to objects with halo masses of 
5-10 X 10"Mq, ages of about 40-50 Myr (cf. Fig.[B, metallicities around 1/3Z0 [cf. eq. (|A9ll] and E(B-V) « 0.15 (cf. 
Fig.|2]l; the predicted mass in stars amounts to 4-5 x 10'' Mq. 

In conclusion, the dust prescription we derive from observations of LBGs is consistent also with the observations of LAEs at 
high redshift. 

4. DISCUSSION AND CONCLUSIONS 

We have extended to the earliest evolutionary phases of high redshift galaxies the physical model worked out by Granato et 
al. (2004) which proved to be capable of accounting for the wealth of data on the heavily obscured active star formation phase 
probed by (sub)mm surveys and on the subsequent passive evolution of spheroidal galaxies. Lapi et al. (2006) have shown that 
this model, which features a coevolution of galaxies and active nuclei, can also provide good fits of the cosmic epoch dependent 
optical and X-ray LFs of active galactic nuclei. 

We have directly borrowed from the Granato et al. (2004) model three basic ingredients: (i) the formation rate of galactic 
halos; (ii) the star formation history as a function of halo mass and redshift, and the IMF of stars; (iii) the chemical evolution. 
These ingredients allow us to obtain the unattenuated spectral energy distribution as a function of galactic age. For the galaxies 
of interest here, which just begin to form their metals, we need a specific recipe for the attenuation. Based on simple physical 
arguments, we have adopted a power-law relation between UV attenuation, SFR and metallicity [eq. (O], tuned to fit the UV 
luminosity vs. E{B-V) correlation found by Shapley et al. (2001); this UV attenuation turns out to be increasing with galactic 
age and with increasing galaxy mass. With this additional ingredient the model accounts also for the observed LFs of LBGs at 
different redshifts. The larger attenuation for more massive objects implies that the UV LF evolves between z « 6 and z « 3 
much less than the halo mass function. This is due to the fact that the higher dust attenuation mitigates the fast increase with 
decreasing redshift of the number density of massive halos. As a consequence, the UV LF is only weakly sensitive to the values 
of cosmological parameters, and in particular of erg. However, if LF estimates will be extended to z ~ 10, the effect of erg could 
be visible down to small enough masses for which the effect of attenuation is small. 

The derived UV LFs imply that LBG can account for a complete reionization at z ~ 7. Most ionizing photons are produced by 
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Fig. 8.— Left: Cumulative distribution of SFRs for LBGs with -24 < M1350 < -18 (dashed line), LAEs with lO^^^ erg s"' < Lljc < 3 x lO''^ erg s"' (dotted 
line), and submm bright sources with 5'85[i,in,^(i+-) > 1 mjy (dot-dashed line) at different redshifts. 

galaxies with halo masses in the range lO'^^M© < Mh < 10" '*Mo at z < 8, that dominate the faint portion of the UV LF. 

We expect that LBGs with detectable Lya emission are on average much younger, less massive and less dusty, in agreement 
with the analysis of z ^ 4 LBGs observed in the GOODS-S survey by Pentericci et al. (2007). Note, however, that a second Lya 
bright phase may be expected at much later times, when the interstellar medium is being swept away by super-winds driven by 
the quasar feedback, allowing Lya photons to escape from the galaxy (see also Thommes & Meisenheimer 2005). 

Only a small fraction (2-4%) of LBGs at z « 3 are detected by deep X-ray surveys (Laird et al. 2006; Lehmer et al. 2005). 
X-ray luminosities estimated using stacking techniques are generally below 10"*^ ergs"', consistent with the luminosities due to 
the starburst activity (Type-II SN remnants, hot gas, and high mass X-ray binaries). The small fraction of z « 3 LBGs with 
Lx > 10^^ ergs"' generally shows properties of active nuclei. A single faint active nucleus was found by Douglas et al. (2007) in 
a large sample of z > 5 LBGs. Our model implies that intense star formation activity is accompanied by an exponential growth 
of a supermassive BH at the galaxy center The corresponding nuclear X-ray emission (cf. Granato et al. 2006) can be written as: 

T Q in44 V fKxV^ ( M, \ _i 

Ly«3x10 — ^ — r erg s , (12) 

^ 0.1 \20) \Mq yr"' ) ^ 

where -q is the mass-to-radiation conversion efficiency, Kx is the bolometric correction for the 2-10 keV band and M, is the 
accretion rate onto the supermassive BH. In the LBG phase the nuclear BHs are still very far from their final masses (proportional 
to the halo masses) and only those associated to halos at the upper limit of the range represented in the available samples 
{Mii > IO'^Mq, see Fig. |3j can reach Eddington limited accretion rates large enough (M, > 0.003 Mq yr"') to yield Lx > 10"*^ 
erg s"' (see Fig. 1 of Lapi et al. 2006); the SFRs of these objects exceed lOOMoyr"' (Fig. [Til. Interestingly, the model predicts 
that in more massive primeval galaxies during the late SMB phase, the SFRs may reach thousands Mq yr"' and correspondingly 
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the nuclear accretion rates may attain tens M© yr~'. In fact, Granato et al. (2006) have shown that the model successfully 
reproduces the X-ray observations of SMB galaxies by Alexander et al. (2005a,b) and Borys et al. (2005). 

We find that the model reproduces the Lya LFs at z = 4.9, z = 5.7 and z = 6.5 (Ouchi et al. 2003; Shimasaku et al. 2006; 
Kashikawa et al. 2006b). The effect of increasing opacity with increasing mass is even stronger for LAEs and is responsible for 
the weaker evolution of the LAE LF, compared to the LBG's, between z = 6 and z = 3, clearly discernible in Fig.|7] Evidences of 
only weak evolution of the LAE LF between z = 6.6 and z = 3.4 have been noted by Taniguchi et al. (2005; see also Murayama 
et al. 2007). 

According to this model, the galaxies currently selected in Lya (Llyq si 2 x lO'*^ erg s~') are very young (cf. Fig.|7]), with 
typical ages up to few x 10^ yr. The Lya luminosity range sampled by surveys carried out so far corresponds to a narrow 
distribution of halo masses around ~ 10"Mq, and stellar masses ranging from 5 X 10^ M© to lO'^M© (see Fig. [B- For given 
halo mass, Lya-selected galaxies have lower stellar masses than LBGs and their stars are preferentially formed in the densest 
central regions; they are therefore expected to be more compact than LBGs, as is indeed observed (Dow-Hygelund et al. 2007). 
Galaxies with higher halo masses have higher Lya luminosities but are rare because of the low number density of massive halos 
at high redshifts and of their short lifetime in the LAE phase; they are thus easily missed by the available surveys, covering very 
small areas. These rare objects would however have very interesting properties: high Lya luminosity with very young ages, 
low metallicities and correspondingly high Lya equivalent widths, low masses in stars within large halo masses. On the other 
hand, the lower luminosity portion of the LAE LF is dominated by relatively older objects, more chemically evolved and with 
lower Lya equivalent widths. Therefore the extension to lower luminosities of the LAE LF will depend rather critically on the 
minimum detectable equivalent width. Both Lya-selected galaxies and LBGs are associated to relatively small DM halos and so 
that they are not expected to show strong clustering, consistent with the results by Murayama et al. (2007). 

Fig.[8]shows the cumulative distribution of SFRs for LBGs, LAEs, and SMBs at different redshifts. Although we must be aware 
that there is a significant overlap among these populations, some trends clearly emerge. At high z most of the star formation is 
associated to the LBG phase, but a substantial fraction of the galaxies with the highest star formation rates are predicted to 
be dust-obscured (and therefore SMB). This is because, in the present framework, the more massive galaxies have higher SFRs, 
yielding higher metallicities and higher dust attenuations. At redshift z > 6, however, the space density of SMBs is very low. Only 
later on, an increasing fraction of the cosmic SFR occurs in very dusty galaxies. Since SFRs are strongly correlated with halo 
masses, SMB galaxies are expected to show stronger clustering, consistent with the results by Blain et al. (2004). Observations 
of the LBG angular correlation function at z > 4 (Kashikawa et al. 2006a; Lee et al. 2006; Hildebrandt et al. 2007) find evidences 
of substantial clustering, yet not as strong as indicated by data on SMB galaxies. In particular, the clustering was found to depend 
on luminosity, suggesting a close correlation between halo mass and SFR, in keeping with our expectations. 

To sum up, our model establishes a coherent scenario encompassing the variety of galaxy populations selected by different 
techniques: Lya emitters, Lyman break galaxies, submillimeter bright, and passively evolving galaxies. In this scenario, these 
populations correspond to subsequent steps in the evolutionary sequence of high-z galaxies. 

We warmly thank A. Bressan for enlightening discussions, and the anonymous referee for very useful comments and sugges- 
tions. This work is partially supported from ASl, INAF and MIUR grants. 
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APPENDIX 

A. A SIMPLE RECIPE FOR THE STAR FORMATION RATE IN PROTOGALAXIES 

The history of star formation, of fueling of the mass reservoir around the active nucleus, and of accretion into the central BH 
can be easily computed by numerically solving the set of differential equations written down by Granato et al. (2004; see also 
the Appendix A of Lapi et al. 2006). In this Appendix we summarize the results concerning the SFR, which are the relevant ones 
for this paper, and, whenever possible, we give simple analytical approximations for them. 

At the virialization time a DM halo of mass Mh hosting a galaxy contains a mass Minf(O) = /cosmA^// of hot gas at the virial 
temperature, /cosm ~ 0.18 being the mean cosmic baryon to DM mass-density ratio. The gas cools and flows toward the central 
regions of the halo at a rate 

Mi„f = -M,o„d-MiT ' (Al) 

where Mcond = Mnf/fcond, and the 'condensation' timescale fcond = inax[fcooi(^H), fdyn(^H)] is the maximum between the dynamical 
time and the cooling time at the halo virial radius Rh- The cooling time includes the appropriate cooling function (Sutherland & 
Dopita 1993) and allows for a clumping factor C in the baryonic component. The second term on the right hand side describes the 
influence of the quasar kinetic energy output on the hot gas distributed throughout the DM halo. Note that in the above equation 
the effect of the angular momentum is neglected, since it is lost by dynamical friction through mergers of mass clouds on 
a time scale fop ~ 0.2((^/ln(^)fdyn, where ^ = Mh/Mc (e.g.. Mo & Mao 2004); major mergers, which are very frequent at high 
redshift, imply ^ ^ a few. 

The mass of cold gas is increased by cooling of the hot gas (Mcond), and decreased by star formation (M,t) and by the energy 
feedback from SNae (M^^^j) and quasar activity {M^^^)- Its evolution thus obeys the equation: 

Mcoid = Mcond - ( 1 - 7^)M. - M,^l - Mfjj° , ( A2) 

where TZ is the fraction of gas restituted to the cold component by the evolved stars, amounting to 7?, « 0.3 under the assumption 
of instantaneous recycling. Strictly speaking, this value of TZ is an upper limit, since only a fraction of evolved stars have a 
significant mass loss in the evolutionary phases considered here; the results, however, are only very weakly sensitive to the 
chosen value, in the allowed range. The cold mass ending up in the reservoir around the central super-massive BH is lost at a 
very low rate (see Granato et al. 2004) and has been neglected in the above equation. 
The energy feedback from SNae is parameterized as 



Mfl = PsnM, , (A3) 

where the efficiency of gas removal 

Nsn^snEsn _^^^f NsN \ f <^sN\f EsN \ f Mh V^^^ f l+zY^ 
= " U X 10-3/Me j loos J llO^ j [rn^J [—) ^^'^ 

depends on the number of SNae per unit solar mass of condensed stars Nsn^ on the energy per SN available to remove the cold 
gas esnEsn, and on the specific binding energy of the gas within the DM halo, fibind- Following Zhao et al. (2003) and Mo & 
Mao (2004), the latter quantity has been estimated for z > 1 as ^bind = V^/(c)(l +/cosm)/2 w 5.6 x 10^'^ (Mh/ 10^^ Mq)^'^ [(1 -I- 
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z)/7] cm^ s~^; here Vh is the halo circular velocity at the virial radius and /(c) sa 2/3 + (c/21.5)°^ ^ 1 is a weak function of the 
halo concentration c. We adopt the strength of SN feedback esn = 0.05 used by Lapi et al. (2006) in order to reproduce LFs of 
galaxies and quasars at high redshifts, and also to reproduce the fundamental correlations between local elliptical galaxies and 
dormant BHs. 
The cold gas turns into stars at a rate 

= / — - « -— , (A5) 

J max[fcool,fdyn] U 

where now fcooi and fdyn refer to the mass shell dMcoid, and U is the star formation timescale averaged over the mass distribution. 
Eqs. ( lAll i and iA2\ can be easily solved setting sa A^coid/f* and neglecting the effects of the energy injected in the gas by 
the accretion onto the central BH. With this approximation the infalling mass declines exponentially as Minf(f) = Minf(0)e~' 
while the SFR evolves according to 



M.(0: ^-^'^ 



fcond(7-l/'S) 



-'Ac™d_„-.!7'/'c<: 



(A6) 



with 7=1 - TZ + (3sN- The quantity s = fcond/f* is the ratio between the timescale for the large-scale infall estimated at the virial 
radius, and the star formation timescale in the central region. An isothermal density profile yields i « 5. 

Although eq. ( IA6l l has been obtained neglecting the quasar feedback, it turns out to be, at the redshifts relevant here (z > 2), 
a good approximation of the results of the full numerical integration of the system of differential equations given by Lapi et al. 
(2006), provided that: (i) we take into account that the star formation is strongly suppressed by the quasar feedback after a time 

Aw«2.5xl0« fil^) yr, (A7) 



\ 1 J 

where Tix) = 1 for jc > 1 and Tix) = for x < 1; (ii) we adopt the following expression for the condensation timescale 

In the latter expression, the scaling with the redshift mirrors the dependence of the cooling and/or dynamical time on redshift. 
Moreover, the weak dependence on Mh renders the impact of the energy feedback from the quasar on the infalling gas; its effect 
is mimicked by allowing longer condensation timescales for the gas in more massive halos, exposed to stronger quasar feedback 
from more massive BHs. 

Note that the SFR evolves differently for different halo masses, see Fig.[T] For large masses Mh 1O'^M0, it increases almost 
linearly with galaxy age in the initial stages t <C fcond, and then it is suddendly stopped by the energy feedback from the quasar 
at f « Afburst. On the contrary, for low-mass halos Mh < lO'^M© it is first almost constant and then slowly declining; indeed, 
in such low mass galaxies the quasar feedback is very mild because the mutual effects of star formation and gas flow toward 
the reservoir around the central BH result in a strong decline of the BH mass with decreasing halo mass [e.g., Shankar et al. 
2006, their eq. (16)]; thus the SFR can proceed for much longer times (note from eq. (A7) that Afbmst ^ 3 Gyr at z w 6 for 
Mh ^ 10" Mq). Note that Fig.[T]refers to redshift z = 6; at z w 3 the SFR retains the same time dependence as at z = 6, but is on 
average lower by a factor of « 2 for a given halo mass. 

The total mass M^""' cycled through stars during the time Afburst can be estimated by integrating the SFR [eq. (IA6I 1I. The 
present-day total mass of surviving stars M™* = /surv^^^"*' can then be derived by assuming a specific IMF; /surv is around 60% 
after 10 Gyr from the burst for a Salpeter IMF, and is around 30% for the IMF adopted here. The latter is a double power-law 
with slope 1.25 from 120Mq to IMq and 0.4 from IMq down to OAMq (Romano et al. 2002; for a review on the IMF see 
Chabrier 2005). 

An interesting outcome of our model is that the fraction of the original baryons associated to the galactic halo that is trapped 
in present-day stars stays almost constant at M"°*/Minf(0) « 0.2 for Mh ^ 3 x IO^Mq, and decreases very rapidly below this 
threshold. As discussed by Shankar et al. (2006), this is in keeping with results obtained through X-ray and weak lensing 
estimates of halos around galaxies. 

When computing the metal content of the cold gas, one has to take into account the enrichment of the primordial infalling gas 
due to earlier generations of stars, and the gas outflows due to the energy injection by SNae and quasars. In the first stages, the 
chemical enrichment of the cold gas component is very rapid; e.g., 1/100 and 1 / 10 the solar abundance in the gaseous component 
is attained after 1.2 x 10^ and ~ 4 x 10^ yr, almost independently of the halo mass and redshift. Only less than 1% and 10% 
of the total stellar mass has been built up at these times and therefore the expected final mass fraction of metal-poor stars is tiny. 
This rapid enrichment is due to the short lifetime of massive stars > IOM0 (that are efficient metal producers) relative to the 
timescale tcon^ of the infall of the diffuse medium with primordial composition, which dilutes the cold star-forming gas. A simple 
approximation of the subsequent evolution reads 

« 10-2 ^ ( ^ ) , for 7 X lO^yr < f < min[Afburst, Af^.t] , (A9) 

7 X 10' yr \ 7 / 

where Afsat = 7 x 10^(Mh/101°Mq)2/3 yr is a characteristic time after which Z keeps almost constant. 
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B. REIONIZATION 

In this Appendix we briefly describe the computation of the reionization redshift predicted by our model. The average global 
production rate of ionizing photons per unit volume at the cosmic time T is 

{Ki){T) = J dNgnN<)nm9n,T) , (Bl) 

where ^912 is given by eq. ^ and its distribution function <l>(A'9i2, T) can be derived in strict analogy to eq. (HJi. The transition 
from a neutral to a reionized IGM can be described in terms of the evolution of the volume filling factor of Hll regions, Qun, 
which is ruled by the equation (Shapiro & Giroux 1987; Madau, Haardt & Rees 1999) 

dQmi ^ (Ngn) Qmi . ^^2) 

dT «H free 

here «h is the comoving HI number density, and the recombination time is given by (Madau, Haardt & Rees 1999) 

1+Z\ I C 



'- = 0-^— j UJ ^^^^ 

in terms of the clumping factor of the IGM, C. The IGM clumping factor has been investigated through hydrodynamical 
simulations (e.g., Gnedin & Ostriker 1997). Iliev et al. (2006) proposed an analytic description of the simulation results: 
C{z) = l7.6e-o i:+oooii; in the relevant redshift interval 7 < z < 10 one finds 10 > C > 7. Here we adopt C = 7 as a reference 
value. Reionization is completed at the redshift when Qhh = 1 (right panel of Fig. |6l). Finally, the optical depth to electron 
scattering is given by 

res(z)= / dz' '^^ 
Jo 



dz' 



cajrieiz'), (B4) 



where rig = Qhii «b (1 + z)^ is the electron density, ng « 2.5 x 10 cm ^ is the present day baryon density, and ctj « 6.65 x 10 
cm^ is the Thomson cross section. 



